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Abstract. In the framework of transformation optics, we show that the 
propagation of a locally superluminal refractive index perturbation (RIP) in a 
Kerr medium can be described, in the eikonal approximation, by means of a 
stationary metric which we prove to be of Gordon type. Under suitable hypotheses 
on the RIP, we obtain a stationary but not static metric, which is characterized 
by an ergosphere and by a peculiar behaviour of the geodesies, which are studied 
numerically, also accounting for material dispersion. Finally the equation to be 
satisfied by an event horizon is also displayed and briefly discussed. 
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1. Introduction 

Effective geometries for light, first introduced by Gordon [1 and extended to 
nonlinear electrodynamics, are able to provide analogue black hole metrics and 
also the possibility to perform experiments involving Hawking analogue radiation 
[l[3ll[Slll[71[Hlini[ini[IIl[Il[IllIll- Recently Philbin et al. proposed an optical 
analogue in which a soliton with intensity /, propagating in an optical fibre, generates 
through the nonlinear Kerr effect a refractive index perturbation (RIP), Sn = 712/, 
where n2 is the Kerr index |15| . The same mechanism has also been generalized to 
a full 4D geometry by Faccio et al. [TB]. The RIP modifies the spacetime geometry 
as seen by co-propagating light rays and, similarly to the acoustic analogy, if the RIP 
is locally superluminal, i.e. if it travels faster than the phase velocity of light in the 
medium, a trapping horizon is formed and Hawking radiation is to be expected. A 
full analysis for the case of a static black hole metric has been recently discussed in 
[TT] . Our aim here is to tackle the following aspects: (i) provide a general derivation 
of the metric in the eikonal approximation from relativistic nonlinear electrodynamics 
in a medium and (ii) perform a first analysis for a generic RIP which gives rise to a 
stationary (but in general non static) metric. The metric we obtain allows to determine 
in a straightforward way the limits under which an ergoregion occurs. Then a study 
of the behaviour of null geodesies follows, with the aim of identifying conditions such 
that trapping of light occurs inside the dielectric perturbation. Finally, the equation 
to be satisfied by an event horizon is provided. 

2. Effective geometry 

In this section, we introduce the specific analogue model that we have analysed. As 
we need to work both in the laboratory frame as well as in an inertial reference 
frame which is moving at relativistic speed with respect to the original frame, it is 
convenient to adopt a covariant formalism. Then, consider a reference frame in which 
the dielectric medium is moving with four- velocity m'^. Moreover, consider a medium 
for which the permittivity and the permeability have respectively the form 

e"^ =£(£;) (7/"'^- u^u-^), (1) 
fi'^^^^ioMv^^-u'^uf'), (2) 

where ry^jy is the Minkowski metric tensor. That is, the permittivity is a function of the 
invariant electric field amplitude E := -y/— i?^£'^ alone. In the laboratory frame, this 
reduces to the assumption that a polarisation is induced in the medium only when it is 
exposed to an electric field, and that the nonlinear effects depend only on the intensity 
of the latter. The permeability is assumed to be constant and equal to the product /i 
of the vacuum permeability /iq the relative permeability fir of the medium. These 
conditions indeed stand true for diverse nonlinear media. For example, for the case of 
an isotropic x'-'^^' medium, in the laboratory frame 

£==eoe. -eo(l + X^''+X^'^^'), (3) 

where Eq is the vacuum permittivity. 

Nonlinear effects are usually very small and a very strong field is needed to 
exploit the nonlinearity of the medium itself. Thus, we assume that there exists 
a strong, background field E, which fixes the properties of the dielectric. Note 
that these assumptions are rather general. In fact, Eqs. (fTl) and (pi) hold for any 
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nonlinear, homogeneous dielectric, no matter how the strong, background pulse is 
generated We are interested in the small and rapidly oscillating fluctuations of 
the electromagnetic field on top of such a backgrounc||] 

According to the results of Refs. [SI [TTl |T2], these small fluctuations, "feel" an 
effective curved spacetime. The result is what a relativist would call a "bi-metric" 
theory. In other words, there exist two classes of rays, namely ordinary and extra- 
ordinary rays, which correspond to two different polarizations of the background 
field. This phenomenon is analogous to birefringence in uni-axial crystals, where 
the ordinary and extraordinary rays each obey distinct quadratic dispersion relations. 
In the analogue gravity terminology, each ray "feels" a different geometry, or analogue 
metric: the Gordon's metric 

1 " 



for the ordinary rays, and the metric 



1 



(4) 



1 - 



1 



7i2(l+$) 



(5) 



for the extra-ordinary rays. In (|4])-([5|, we have introduced the "nonlinear parameter" 
^ := the background field polarization four- vector and the refractive 

index n := ^er[ir- In cartesian coordinates (t, x, y, z), 77^1, = diag(c^, —1, —1, —1). 

In [Appendix A| we show that, under reasonable approximations, the extra- 
ordinary metric is equivalent to the Gordon one, albeit with a different effective Kerr 
perturbation of the refractive index, i.e. n = no + 3o?n instead of n = no + dn. Thus, 
we restrict our attention to the Gordon metric. 

In the following we study the lightlike geodesies of the analogue metrics. It 
is important to note that these geodesies describe the spacetime trajectories of 
monochromatic waves. In the present paper we often refer to these trajectories 
with the term "light rays". However these should not be confused with the 
usual light rays referred to in geometrical optics which, in the Newtonian picture, 
describe the trajectories of wave packets (or "photons"). For this reason, in the 
forthcoming sections, we analyse the monochromatic components (congruences of 
lightlike geodesies) of the electromagnetic field and the analogue geodesic motion of 
the related constant phase surfaces, which are precisely the objects involved in the 
study of the analogue Hawking radiation [18l [19]. By contrast we do not consider the 
notion of "photon" (or wave packet) here. 



3. The metric in the RIP frame 



Consider the case of a pulse moving with uniform velocity v along the positive x axis. 
According to the discussion above, the pulse excites a localized RIP which also moves 
accordingly. The latter can then be analytically described as a functio n of cc — vt and, 
assuming cylindrical symmetry, of the transversal coordinate p = yj-ip + z'^ . Thus, 
in the laboratory frame, n = no -I- 5n{x — vt,p). In this situation, it is convenient to 
express the Gordon metric in a reference frame which is co-moving with the RIP itself: 

ds^ ^ 0^72 1 + dt'^ + 272^ (1 - n2) dt'dx' 

n^V c/V c / n^ 

f For reasons connected with the actual experimental verification of this theoretical framework, we 
may think of this strong background field as a spatially localized (albeit rather extended) pulse, which 
is travelling through the medium. 
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-j^(l + -)(l--)dx'^-dp'-p'dcl,^ (6) 
V nc/ V nc/ 

where the primed coordinates are relative to the co-moving frame, and 

n := riQ + dn{x' , p) . (7) 

Here, Sn is a C°° function, rapidly decaying at infinity and with a single maximum of 

height rj, describing the RIP. It is then evident that both dt' and 9^ are Killing vectors 

for the given class of metrics. Any specific choice for the function Sn gives rise to a 

specific metric in the aforementioned class. We point out that an isotropic refractive 

index in the laboratory frame corresponds to an anisotropic refractive index in the 

pulse frame due to length contraction associated with a boost in the x direction. 

If the refractive index depends explicitly on p, the metric is stationary but not 

static: the integrability conditions stated in Frobenius theorem [20] are not satisfied. 

In particular, if = (gooi 5oij 0; 0) are the covariant components of the timelike 

Killing vector dt', the integrability conditions £,[x'^fi^iy] = 0, for A = 0, /i = 1 give 

?o(ViCi/ - Vj^Ci) + ^i(V,,Co - VoCi.) = 0, i.e., for i' = p, 5oo(-9p.goi) + goidpgoo = 0, 

which is easily seen to be not satisfied. The static case is obtained if Sn does not 

depend on p, and it has been extensively studied in [17 . Here, we consider the more 

general stationary, non-static casej|] We first point out that for our class of metrics 

an ergosurface is allowed, given by the locus of points satisfying the equation 

nv , , 

500 = ^ 1 =0. 8 

c 

This equation admits solutions only for the range of RIP velocities v given by 

\ V \ 

< - < — . 9 

no +r] c no 

At the ergosurface, which is the boundary of the ergoregion, the Killing vector df 
vanishes and upon entering the ergoregion it changes sign. Here, in contrast with 
the static case, the ergosurface does not coincide with the event horizon. Indeed, the 
given hypersurface is not a lightlike one, as can be easily verified by looking at the 
behaviour of its normal vector. However, for p = it actually becomes lightlike, i.e. 
it meets the event horizon in the antipodal points which solve 

Sn{x\0)^^^^. (10) 

V 

For stationary nonstatic situations it is more difficult to identify a nontrivial event 
horizon in effective geometries (see e.g. Visser's analysis for the case of acoustic black 
holes [H]). Before discussing the question of the existence of an event horizon, we 
investigate whether the path of light rays, that is the null geodesic curves of metric 
do suggest the existence of such a horizon. 

§ A simple, non-technical definition of static and stationary spacetimes can be given as follows. A 
spacetime is stationary if it is possible to find a reference frame in which the metric coefficients do not 
depend on time. This means that time translation is a symmetry of the spacetime. However, such 
a spacetime might not be invariant under time-reversal. In the case of time-reversal invariance, the 
spacetime is not only stationary, but also static. An example of static spacetime is the Schwarzschild 
solution, describing the static empty spacetime outside a non-rotating black hole. If the black hole has 
a non-vanishing angular momentum, the spacetime is still invariant under time translations, however, 
time reversal reverses the sign of the angular momentum, turning a clockwise spinning black hole to 
an anticlockwise one. Then, the spacetime around a spinning black hole is stationary, but not static. 
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4. Geodesies and trapping of light 

In the laboratory frame the metric (|6| is 



-dt^ - dx^ - dy^ - dz^ 



(11) 



■n?{x — vt, y, z) 

We model the RIP with a gaussian shaped function which propagates along the x axis 
with velocity v. We also suppose that, during the propagation, it keeps its shape both 
in space and in time. Assuming cylindrical symmetry, we need to account only for 
one of the two transverse coordinates, y, say. Then, we write the refractive index in 
the form 



no 



(5n, 



5n — rje 



Thus, the geodesic equations for null rays take the explicit form 
2 Sn 



-{x - vt)cH^ 



2 6r 
Sn 

n 2 



rye t 



v{x — vt)P — yyi — {x — vt)xi 



(12) 

(13) 
(14) 
(15) 



where the dot denotes derivation w.r.t. an afline parameter A. A trivial solution is 
^ = c/n, y(t) = 0. However, in the general case a numerical analysis is required. 
In the following subsection we solve the equations above using the Matlab's ode45 
function. This analysis will provide an insight in the underlying physics. 

4-1 ■ Numerical solutions and trapping of light 
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Figure 1. (a) Sketch showing the initial conditions of the simulation, (b) Light 
ray phase velocity (c/n). The shaded region highlights the range of velocities of 
the RIP for which the rays are trapped. The two straight lines indicate the values 
of v/c used in the simulations. 



Along a geodesic we have {ci,x,y,z) = {k^ , , k"^ , k^) = {n'^oj/c,kx,ky,kz), and 
taking the first derivative (with respect to the affine parameter A) of the relation 
ct = n^w/c, we find an equation for uj which determines the evolution of the frequency 
of the ray along a geodesic. At this point a comment is in order: here we have 
defined uj as the zero-th covariant component of the four-momentum k^. It is easy 
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to sec that this corresponds to measuring frequencies as the inverse of the lapse of 
coordinate time between two successive crests (this time is the Minkowskian time 
measured by the laboratory observer). This is of course different from the analogue 
proper time dr — ds/c of the analogue metric ( |11[ ), which has no physical significance: 
experimental measurements refer to the true Minkowski metric and not to the analogue 
curved metric. 

Then, for the evolution of the frequency u along the path of the ray, we obtain 
t — 2(jjn5n 

uj= . (16) 

For the sake of definiteness assume the propagation to be in fused silica and the rays 
to have initial wavelength Ai„ ~ 527 nm. We assume the RIP to have cr = 1 m and 
amplitude rj = 0.1, which is a quite unrealistic situation, but simple to deal with. 
Indeed, the aim of the simulation is merely to study the qualitative behaviour of the 
geodesies. We have verified that this behaviour does not change even for much smaller 
and realistic values of cr ^ 1 — 100 fim. 

The other parameters of the simulation are the velocity v of the RIP and the initial 
conditions of the rays, i.e. the array {Un, Ximyimtin, iimVin}- In the simulation, 
the rays start away from the RIP (which is right moving) and their initial velocity is 
directed along the x axis, as shown in the sketch of Fig. [T| a). More precisely, the initial 
conditions are {0, y™, ti„, ii„, 0}, with the additional constraints iin/Un = c/no 
and din = "o'^m- 

The results of the simulation show that the behaviour of null geodesies is strongly 
dependent on the velocity v of the RIP. Namely, there exists a range of velocities for 
which the rays are trapped [ji] and these coincide with those predicted by Eq. ^ 
(shaded region in Fig. [l|b)). On the other hand, for lower values of v, the rays 
which enter the RIP are bent and eventually escape. In Fig. [l|b), the black curve 
represents the velocity of an axial ray, i.e. of a ray propagating along the x axis, 
versus x — vt. In the same figure, we represent with horizontal lines two choices for 
the RIP velocity, representing two possible situations: a faster (green line) one and a 
slower (red line) one. In the first case (faster RIP), the green line intersects the black 
curve. As the ray approaches the point corresponding to the first intersection it slows 
down due to the increasing refractive index, until it acquires the same velocity of the 
RIP. Asymptotically, the ray is "trapped", and moves together with the RIP. In the 
other case (slower RIP), the ray velocity is always larger then the RIP's one (red line 
of Fig.[l];b)). Then, the ray can cross the RIP and escape. In this "subluminal" case, 
no trapping occurs. 

Subluminal RIP The RIP propagates at v — 0.56c (red line in Fig. [l|b)), and it is 
always slower than the light rays. In Fig. [2][a) we show the numerical results of the 
simulation in the co-moving reference frame. The rays start on the left of the RIP. 
They propagate along straight lines with velocity c/uq until they reach the RIP. Then, 
inside the RIP, they slow down, bend and finally escape from the RIP, propagating 
along straight lines in different directions, depending on their impact parameter 
The projection of the trajectory on the {t',x') plane, which is also shown Fig. [2]ja), 
evidently shows a change of slope corresponding to the trajectory being inside the RIP. 
This means that the rays, while crossing the RIP, slow down. Fig.[2]jb) represents the 
evolution of cj along the x direction in the laboratory reference frame. We observe 

II This is true also when taking into account the material dispersion, which is neglected in this 
subsection. 
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Figure 2. Numerical results in the co-moving reference frame (a) and frequency 
evolution in the laboratory reference frame (b) for the subluminal case. 

that, as expected, the frequency varies consistently when the rays reach the RIP. The 
rays go through a blueshift inside the RIP and then they emerge with a Doppler 
redshift, due to the fact that the rays are bent, and thus emerge at different angles. 
This behavior is described in more detail in [IB] . 

Trapping RIP: The RIP propagates at v = 0.66c (green line in Fig. [ijjb)). In 



xio" 




Figure 3. Numerical results in the co-moving reference frame (a) and frequency 
evolution in the laboratory reference frame (b) for the trapping case. 

Fig. |3]^a) we show the numerical results in the co-moving reference frame. The 
rays propagate initially along straight lines with velocity c/tiq (for clarity, we recall 
that rays propagate along geodesies with the phase velocity of the corresponding 
electromagnetic field) . When they reach the RIP, they are bent and at the same time 
slow down. As shown in Fig.[ljb), at a certain point, away from the center of the RIP, 
the rays will have exactly the same velocity as the RIP, and can no longer escape from 
it. In the co- moving reference frame they continuously slow down and at the same 
time they oscillate in the transverse direction with decreasing amplitude oscillations, 
until, asymptotically, all of them stop at the same point. This point, represented in 
red in Fig.[l]Ja), is the white hole horizon X- (see next section). If we look at the slope 
of the projection t' — t'{x') we see that it diverges, which means that the rays velocity 
tends to zero corresponding to this point. Clearly, in the laboratory reference frame 
the rays do not stop but they move together with the RIP. As shown in Fig. [sf^b), 
as the rays reach the RIP, their frequency starts growing indefinitely, so we have an 
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infinite blueshift. 



4-2. Effects of dispersion 

In the above considerations we did not account for the material dispersion, that is 
we assumed n{T,t) independent of u. However, the refractive index of a medium is 
in general a complicated function of the frequency uj of the incident wave, that is 
n = n{r,t,uj). We restrict our attention to the case of "normal dispersion", that 
is we assume that the refraction index grows with frequency. In this case, the full 
dependence of the refractive index from uj is phenomenologically described by the 
Sellmeier equation. A standard way of treating dispersion is to expand n{uj) in power 
series of uj (about an initial frequency wq), or equivalently, by the relation ck = tun, 
to expand k in power series of uj as follows: 

The material dispersion also affects the metric, which therefore depends on the 
frequency: 

ds^ = — -dt"^ - dr^. (18) 

Therefore, different rays propagating in the medium "feel" different metrics, depending 
on their frequency. Moreover, even a single ray will see an evolving metric while 
travelling, since its frequency (and hence the effective metric) changes during the 
propagation. The geodesic equations ( 13p5 1 keep the same form as in the non- 



dispersive case, while Equation ( |16[ ) now takes the form 

c^t — 2u)n6n ,^ ^. 

n[n + ZLudnQ/duj) 

We highlight that, contrary to the non-dispersive case, the evolution of the geodesies 
depends on the change in frequency of the ray, owing to the explicit dependence on 
UJ in the geodesic equations. In the case of a subluminal RIP, dispersion does not 
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Figure 4. Numerical results in the co-moving reference frame (a) and frequency 
evolution in the laboratory reference frame (b) for the trapping case when we 
account for dispersion. 



modify the qualitative behaviour of the light ray, which remains very similar to the 
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dispersion-less case shown in Fig. [2j On the other hand, for a trapping RIP, the pres- 
ence of dispersion modifies substantiaUy the geodesies. The results of the numerical 
simulation are shown in Fig. [4] The trapping feature of the horizon is still evident. 
However, taking dispersion into account, the rays are no longer indefinitely trapped 
inside the RIP, and eventually escape. This behaviour is clearly depicted in Fig.jijja), 
where the results are shown in the co-moving reference frame. Note that, at the be- 
ginning the (phase) velocity of the rays is positive, which means that the rays are 
faster than the RIP, but then it becomes negative, implying that the rays have be- 
come slower than the RIP. Let us clarify this aspect. The phase velocity of a ray is 
defined as Vr = to/k. Under conditions of normal dispersion, i.e. > 0, as in our 
case, k grows faster than linearly with u (while without dispersion it would have a 
linear growth). When the rays approach the RIP, their frequency grows and, at the 
same time, their phase velocity decreases. Thus, the rays will begin to lag behind with 
respect to the RIP. When they are out of the RIP they will keep their final frequency. 
In Fig.[4]jb) we show that the frequency grows only by a finite amount. Therefore, the 
rays leave the RIP with a finite (albeit very large) blueshift. 

Trapping time By increasing the velocity of the RIP from v — 0.5c up to w = 0.65c 




Figure 5. Pictures (a) to (g) represent the trapping dynamics in the absence 
of dispersion for different RIP velocities: (a) v = 0.5; (b) v = 0.56; (c) v = 0.6; 
(d) V = 0.62; (e) v = 0.63; (f) v = 0.64; (g) v = 0.65. Picture (h) represents 
the interaction time versus the velocity of the RIP in case of no dispersion (blue 
curve) and in case of dispersion (red curve) . 



(i.e. approaching the trapping case), the rays spend more and more time inside the 
RIP. In Fig. [sf^a-g) we show this behaviour for different velocities, in the dispersion- 
less case. In Fig. [5jh) , we have plotted the interaction time versus the velocity of the 
RIP in the dispersion (red curve) and dispersion-less (blue curve) case. Neglecting 
dispersion, the interaction time tends to infinity as the velocity of the RIP approaches 
the trapping case. We can interpret this divergence as a (meta)stable state where the 
photon is trapped inside the pulse. An analytical estimate for this trapping time is 
given, for the axial geodesic by the formula ttrap = | /J^ i-n{i)v/c ' where r represents 
the longitudinal extension of the RIP. In the dispersion-less case ttrap diverges whereas, 
in the dispersive case, it is always finite. 
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5. Horizon in optics 



In this section, we provide a definition for the horizons from the optical point of view, 
in the dispersion-less case. For the sake of definiteness we focus on the "white hole" 
horizon, but similar considerations hold true for the "black hole" horizon. We work 
in the laboratory frame, where the pulse is moving in the x direction, with velocity v. 



The situation is depicted in Fig. 6(a) 



Recall that the ergosphere is defined as the set of points where the velocity of 
light is the same as the velocity of the pulse, which means n{X) = c/v for all points 
X on the ergosphere. In particular, it follows that a light ray reaching the ergosphere, 
for example in P, will not enter the ball if it is moving in the x direction. However, 
if it is directed along a NW-SE direction (for example toward the center) then it 
will cross over the ergosphere thus entering the ball. This is because it "sees" the 
surface move slower in such a direction: indeed, note that a plane forming an angle 6 
w.r.t. X and moving with velocity v in the x direction, appears to move with velocity 
u — wsin^ < c/n{P) along its normal direction, see Fig. |6(b)"| 





(a) Ergosphere and horizon. 



(b) Surface speed 



Figure 6. (a) Picture of the ergosphere and of the hypothetical horizon at a fixed 
point of time in the laboratory frame. The small segments at P and Q represent 
infinitesimal portions of null geodesies (light rays) passing, respectively, through 
P and Q at that given point of time, (b) Normal velocity u compared to the pulse 
velocity, v. 



Therefore, only a single point on the ergosphere, namely the axial point 
behaves as a true white hole horizon. In order to determine whether a full event 
horizon extends beyond we impose the condition u = c/n{Q) on some point 
internal to the ergosphere (where n{Q) > n{P) — c/v) which then gives 



vsine{Q) = c/n{Q). 

Now, suppose that the horizon is described by a surface p — h{x 
radial coordinate. Then 

tan^?(Q) = f'(Q). 
ax 



(20) 

wi), where p is the 
(21) 



Note that Eqs. (201 and (21) are indeed satisfied by the axial point x-^ (for which Eq. 



(21 1 gives 



90°). A simple manipulation gives 
dh c 
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1 
X (arb. units) 



Figure 7. Solid lines: real part of the solutions to Eq. ( |23| . Dashed line: 
the ergosphere predicted by Eq. |9]l. The parameters are no = 1.45, rj = 0.1, 
V = c/1.52. 



Finally, in the comoving frame, using coordinates the horizon equation is 



dh 
dx' 



± 



(23) 



■n^ — c 

The same equation can be obtained in the following way. It is a common strategy 
to parametrise the horizon by means of a level-set function F such that _F = on 
the horizon; the lightlike condition for the event horizon is g^'^ {d^F){di,F) = 0. In 
our case, there are two further requirements: axial symmetry and that the horizon is 
stationary (i.e. dt'F=0). Then, in the RIP frame the condition which imposes the 
normal to the event horizon to be lightlike is: 

„,2 ' 



7' 1 - n\x', p)^ {d,,Ff + [dpFf = 0; 



(24) 



in the static case. 



dpF 



0, and then, from (24) we deduce that the horizon consists 
le X axis, located at the solutions x'^ 



of two infinite extended planes, orthogonal to t 

of 1 — n^(a;')^ = 0, as described in [T7j. In the nonstatic case, it is evident that 
the horizon has to be contained within the ergosurface, because it is necessary that 
1 — n^(a;, p^v^ jt? < 0. Moreover (at least locally, due to the implicit function theorem), 
we assume that F{x',p) = <^ p - h{x') = 0, i.e. F{x',p) = p - h{x'). Then p4 ) 

±l/-\/— which 



becomes g^H^)^ 
explicitly corresponds to (l23|) 



0. As a consequence, we obtain 



dh 
dx' 



According to Eqs. (20) and (21), as described above, the points x_ and x+ for 



p = belong to the event horizons, which then represent initial points for the ordinary 
differential equation ( 23 1 . Allowed solutions should start from the aforementioned 



points and remain inside the ergosurface. According to numerical evaluations, carried 
out for 77 — 1 — 10~^, in the range of interest for our assumption that the RIP represents 
a small perturbation with respect to the refractive index value in absence of the RIP, 
the RIP exhibits a somewhat odd behaviour in the sense that the solution does not 
extend beyond X- and x+. 

Physically, a light ray passing through a point laying outside the ergoregion 
can travel in either direction (toward the RIP, or away from it). Such a ray is not 
"trapped" , and therefore such a point cannot belong to an event horizon. 

Figure [7] shows the solutions to Eq. ( pS]) for the case 77 = 0.1. The dashed 
line shows the ergosurface predicted by Eq. i9l and the solid lines are solutions to 
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Eq. (23 1 for different initial conditions. These curves represent lightlike surfaces 
passing through points inside the ergosphere (in the RIP frame). The different null 
hypersurfaces which pass near a;+ (or X-), cross the ergosphere at a point which is 
progressively closer to x+ (or x-). This is a numerical evidence of the fact that an 
hypothetical event horizon cannot extend beyond the points x+ and X- . As can be 
seen, the solutions that pass through and x_ are tangent to the ergosurface and 
always remain outside the latter so that the horizon effectively reduces to a point. We 
note that strictly speaking even though these points behave like points of a trapping 
horizon, as shown for example in Fig. [3j it is of course not possible to define them as 
event horizons in the sense of a lightlike hypersurface. One may wonder whether a 
"trapping point" is expected to still emit Hawking-like radiation and, if so, with which 
features. In particular, can we rely on some arguments in order to support the idea 
that there is a pair-creation process which, moreover, still shares common features with 
the Hawking effect? We show in a forthcoming work that the answer is affirmative 
[T5] . Points where the ergosurface becomes lightlike, beyond their properties to be 
trapping points, which behave as standard horizons, share the further property to 
involve a production of particles with thermal spectrum, again like in the case of 
standard (non-degenerate) horizons. 



6. Conclusions 



We have derived the effective metric associated with a RIP in a nonlinear Kerr 
medium. Such a metric simulates the effect on light rays of the presence of a dielectric 
perturbation in a sample of an isotropic x^^^ dielectric material. Then, we have 
investigated the possibility for an ergoregion to arise, and the behaviour of geodesies 
by means of numerical simulations, also taking into account the presence of dispersion. 
Finally, we have provided an equation for the event horizon that highlights a rather 
peculiar geometry, as it appears that the solutions may reduce to a single point. 
Numerical evaluation in the range of interest confirms the former conclusion. Still, 
the two special points X- and x+ of the ergosurface, exhibit trapping properties like 
those of a standard trapping horizon from a classical point of view. 



Appendix A. A remark on the extraordinary metric 

We show that, under reasonable approximations, the extraordinary metric can be 
recast in the Gordon form. We work in a reference frame at rest w.r.t. the medium 
(the laboratory frame). Consider a background field E which, in the laboratory frame, 
has polarization vector orthogonal to its direction of propagation x. We can choose 
Cartesian coordinates {t,x,y,z) so that = (0, — 0,ly — QJz = !)• The line 
element for the extraordinary rays takes the form 

dsf. = "^df - dx'^ - dy'^ - ( 1 dz'^. (A.l) 

^ ' n'' V 1 + 4 / 

Note that the factor in front of dz^ cannot vanish. Moreover, ^ <C 1 (see also |T7]). 
Then, we can safely keep only the leading order term in the z part of the metric so 
that, in so doing, the extraordinary metric reduces to 

dsf^) ^ ^dt^ - dx^ - dy"^ - dz^ . (A.2) 
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Under this approximation, the effective metric for extraordinary rays has again the 
same form as the Gordon one, with an effective refractive index 



to be compared to the refractive index for the ordinary rays, which is 

n = VMi-Er = "0 + -; t^ttt = + on. (A. 4) 

I (1 +;^(l))2 

Here, no is the refractive index of the medium when we disregard nonhnear effects, 
i.e. 



Then, the extraordinary rays "feel" a disturbance of the refractive index three times 
stronger than the ordinary ones. The same conclusion can be also achieved in the con- 
text of Hawking radiation production by a RIP [17] . Indeed, it can be shown that the 
differential equation which defines the extraordinary field modes is well approximated 
by the same equation as for the ordinary ones, apart from an effective RIP three times 
larger. 
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